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ABSTRACT 


A proposed neutron-neutron scattering experiment in which the neutron 
source would be colliding beams from an underground nuclear explosion,was 
investigated. In particular, the accuracy to be expected was determined 
so as to allow evaluation of the worth of the experiment. Phase shift 
analysis was performed for Gaussian, Exponential, Yukawa, and Square 
Well spherically symmetric potentials. The energy range considered was 
20 kev to 2.6 MeV. Scattering length was found to be determined to 
within 6% and effective range within 40% for a 10% scattering cross section 
accuracy. Potential depth would be determined to at least .5% and range to 


‚3%. 
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1. Introduction 


A low energy neutron-neutron scattering experiment utilizing а 
single underground nuclear explosion has been proposed and studied. [2,3] 
Preliminary studies indicate the scattering cross section might be meas- 
ured to an accuracy of 10% from 20 kev to 2MeV. The efforts of this 
paper are directed to a numerical investigation of the expected scatter- 
ing with special attention given to: (a) the possible existence of the 
dineutron; (b) the effect of changes in effective range and scattering 
length upon the range and depth of the spherically symmetric Gaussian, 
Exponential, Yukawa, and Square Well Potentials; (c) the effect of changes 
in effective range and scattering length upon the scattering cross 
section from each of the potential shapes in (b); and (d) a comparison 
of calculated cross sections from the four potential shapes in (b) with 
the shape-independent cross section. 

Precise determination of the neutron-neutron scattering parameters 
would be important for a number of reasons as stated by Moravcsik. [3] 
These reasons include: (1) an investigation of charge independence; (2) 
an Investigation of charge symmetry; and (3) the existence or non-existence 
of a slightly unbound dineutron resonance state. 

This investigation attempted to locate a dineutron state in section 
2 to provide a theoretical justification for its experimental quest. In 
Sections 3 and 4, 2 to 16% deviations in the ratio of effective range to 
scattering length from the accepted n-p ratio were investigated for their 
effect on scattering cross sections and Gaussian, Exponential, Yukawa, 
and Square Well Potential parameters. These results should allow some 
determination of the accuracy in experimental results required to 


determine to a given accuracy the shape-independent scattering parameters. 





In Section 5, the cross section from effective-range theory is compared 
with cross sections from potential scattering to test the validity of 
the shape independent approximation. This report concludes with a 
summary in section 6 and the complete computer program used in the 
calculations along with the calculation results as a part of the appendix. 

The analysis was performed using phase shift scattering methods. 
The well established scattering length and effective range for neutron- 
proton scattering, ак. and ЯҒЫ were taken as a base for calculations. 
The ratio of effective range to scattering length was formed to yield 
a spherically symmetric potential range and depth for square well, 
exponential, gaussian, and yukawa potentials. The absolute value of 
the scattering length was then varied in 2% steps up to 16% while 
holding the effective range constant at its accepted value. 17 sets 
of potential parameters were thus obtained for each of the four poten- 
tial shapes. Then, using the maximum and minimum acceptable value of 
effective range as a base, the effective range was varied in 2% incre- 
ments up to 6% while holding the scattering length constant at its 
accepted value. Eight additional sets of potential parameters were 
obtained for each of the potential shapes. 

The spherically symmetric radial wave equation was then numerically 
integrated using a fourth order runge kutta method for all calculated 
sets of parameters and potential shapes. Integration was performed for 
30 values of center of mass energy ranging from 20 kev to 2.6MeV and 
was carried out to a radius at which the effect of the nuclear potential 
was no longer felt. 

The logrithmic derivative of the numerically integrated solution 
was then compared to the logrithmic derivative of the exact field-free 


asymptotic solution outside the range of the nuclear potential and a 








phase shift obtained. 

The total cross section as a function of center of mass energy was 
then calculated from the phase shift for each case. A ratio of each 
cross section to that obtained from the accepted neutron-proton scatt- 
ering potential parameters was formed. This gave the per cent change 
in total cross section for each incremental change in scattering length 
and effective range. 

Per cent change in total cross section from that for n-p scattering 
parameter calculations was plotted against energy for significant incre- 
ments in two different energy ranges, 20 kev to 100 kev and 100 kev to 
2.6MeV. Potential depth and range were plotted against the ratio of 
effective range to scattering length for each potential shape. A com- 
parison of shape-independent effective range scattering calculations 
with potential scattering calculations was made by plotting the total 


cross section for each case on the same graph. 





2. Possible dineutron state 


An attempt was made to justify the possible existence of a dineutron 
resonance state near zero energy using a compound nucleus model, but 
failure was met in attempting to formulate a compound nucleus model for 
neutron-neutron scattering. 

Numerical calculations for neutron-neutron potential scattering 
were made for a large number of potential depths on the Control Data 
Corporation 1604 digital computer in an attempt to find possible 
resonances, particularly a sharp resonance. Gaussian, exponential, 
and yukawa potentials were used with potential depths ranging from 
0 to the value obtained from the accepted n-p scattering parameters 
in 2.5 MeV steps over two center of mass energy ranges, 20 to 100 kev 
and .| to 2.6 MeV. Results were scanned for resonances. No reson- 


ances were expected and none were found. 





3. Effective range and potential parameters 


Detailed comparisons with experiments have given for the singlet 


effective range and scattering length for neutron-proton scattering: 


раз 2.4 64 0.12 ΕΜ о =-23.677 +0.029 FM CA 


For an examination of the effect of deviations from this scattering 
length on the scattering cross section, the ratio ra, /a., was formed 
and a, was varied in magnitude from 0-16% in 2% steps while holding r,,ç 
at Its accepted value, 2.46 Fm. To examine the effect of deviations in 
effective range on the scattering cross section, the ratio Го) в Mo 
was again formed using the accepted value for a¿, -23.677 Fm, and the 
maximum and minimum values for r,. ,2.58 Fm and 2.34 Fm respectively. 
ros was varied up to 6% in 2% steps in the direction of increasing 
deviation from the accepted value, 2.46 Fm. 

Four spherically symmetric potential shapes were then examined. 
These potentials are: 

Square Well: V = - м Ts Tn 


V = 0 i > ra 


Gaussian Well: 


< 
n 
\ 
< 
о 


ex p EC m 
Yukawa Well: 


V= ХА exp С-га] 


= 
_ Exponential Well: / Τη 


М = — Vo exp Es г/п) 





Blatt and Jackson have arrived at a correspondence between the 
ratio of effective range to scattering length and the potential para- 


meters, Vo and r for each of the well shapes investigated. [1] 


n? 
The well parameters are given in terms of the intrinsic range, b, and 
the well depth parameter, S. The well depth, Vir)» for the potential 
in question is multiplied by a constant such that the first resonance 


occurs at zero energy. The well depth parameter is then the inverse 


of this constant given by: 


Vin = εν (гу 


RES 
The intrinsic range, b, is determined by the shape and range of the 
adjusted potential through the radial wave function for that potential. 


It is defined by: 


= - = Ë. 
ben = г с! DT dr 


RES 
Blatt and Jackson provide interpolation formulas to relate V, and 
ra to rg and a. The form of their formulas have been changed slightly 
to fit the definition of the potentials in this report. іп their 


modified form, these interpolation formulas are: 


Square Well: 
ба С1- 1,574. /а,) +.353 (гъ/03 1 /1- 2.394 0 /ау + 1.80 (ге/о2 | 


b-f/[l-.500(S-I). 486 (s- I= .476(8-1)° πο. Ta(s-l)-. 15% (s-1y 7 


- 24 
Y = 1.022 76x10 s/b? 


In = b 
Gaussian Well: 
$z2Cl-2,030 (4/0) 1.587 (R/O 3 /Cl- 2910 (n/0)4+1.805 Сго/дУ 2 


БЕ гЬ/С!-. 641 съ-1) 1.56& (S-1) — Tel (S-Aj 412 cs--i.eists-i ] 


-24 
V= 2.294208 xı0o s/p? 





Lopes. 06047% 


Exponential Well: 
c 
$zCI- ?.47139(r5/0o) * 1.3 602 (ro/aoS 3 /С | - 3.7134 0 /0.) + 5.1177(г/а) ] 


Бе Г/Г!-.9046-0 +.145 65-19 - 543 [8-19 4.174 (9-1 
-24 
Vo= 7.51541 X10 S/b* 


mz: b/3.5412 
Yukawa Well: 

[4 
S2Cl- 3.5093 (%/0) 42.0 354 (πι) α/{ι- 4.134 6 (го /0-у + 3, 7072tre/0) ) 


ар 
bs m /L1-1.369 ($-1) & 1.043(5-D^- 11e n(s-I #1, ообъ- 1 -,Зг4(- 2 
Vo = 3.1282 1166 χιο τς / εξ 


In = b/2,1196 
The results of these calculations are presented both in tabular and 


in graphical form in tables | and I! and in graphs 1 through 16. The 
tables are arranged to exhibit V9 and га as a function of the percentage 
deviation of га and a from De d and a*., the best values from n-p 
scattering. The graphs depict Vo and га as a function of ro/a using 
the left and bottom scales respectively. Using the top and right 
scales, the percentage deviation of Vo and rp from the best n-p scatter- 
ing parameters is shown as a function of the percentage deviation of 
ro/a from r*, s/aws. 

Variations in potential parameters were found to be of the order 
of a 5% increase in [Vol and a .5% decrease іп г. for a 10% increase 
in scattering length. A 10% increase in effective range, however, was 
observed to cause a decrease in |Volof approximately 8% while increasing 


гп about 9%. This data is shown in table 111. 
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TABLE 111! 
PERCENTAGE CHANGE IN WELL PARAMETERS CAUSED BY A TEN PERCENT INCREASE 


IN SCATTERING LENGTH AND EFFECTIVE RANGE 
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h. Phase shift analysis 


A partial wave analysis of neutron-neutron scattering was performed 
for all potential parameters calculated in section 3. 


Assuming a spherically symmetric potential, the radial wave equa- 


tion 
Та Зе (СВ) На ЕУ - ALP R= 0 (1) 
becomes 
e 
- ὅς SH νο + AT = EU (2) 


where Ц > + г REY 
Rearranging, equation (2) becomes 
ou =-{RrVin- 110) U (3) 


where k^ 2 ME/t? * [EtveD / 41.449] 1079 


апа \/ г) = МУСИ егумен) 4.490102 [63 


Ms neutron-proton reduced mass 
Ес center of mass energy = 1/2 laboratory energy 
Separating equation (3) into two first order equations for simul- 


taneous numerical integration gives 
Z= (4) 
. dr 


— в | 
4 - t V (r> — LED Uu. (5) 
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Equations (4) and (5) were numerically integrated using a fourth 
order Runge-Kutta method with Boundry conditions: 
For ΓΞο 
цо) = О 


Z (o) = 1 


Integration was carried outside the range of the nuclear potentials 
with 500 to 750 steps to about 15 fm for two center of mass energy 
ranges; 0-100 kev in 20 kev intervals, and .1-2.6 MeV in .1 MeV intervals. 
The values of Z and U at r = 15 fm are referred to as Z, and U,. 
The exact asymptotic solution of the field-free wave equation, 


equation (1) with V' = 0, is given by Zeleny as: [77 

Кус =Аз {cos g jn — sin 6, ny (ar)? (6) 

from which 

q = Ag {cos s,[ Bj, - +0) V (Rr) sm s[-&, (an, can (7) 
— (0 ny (kO) J f 


With the substitutions 
0, =r Rp 


dur ord 
dr T +h 
equations (6) and (7) become: 


U= Arr ECOS Se j (kD = SIn 5, п, суў (8) 


240. Ay (tl By cos &, (4 )ل‎ N) μμ (кг) (9) 


А 910 6, (4 omo (КО) 
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where Š, is the scattering phase shift and Jp and n, are spherical 


4 
Bessel and Neumman functions respectively. 

The scattering phase shift, S, , is obtained by matching the 
logrithmic deviations of the exact field-free and numerical solutions 
at a range where the nuclear force is negligible. This procedure 
yields; 

For 2 = 0 
Zn| = £rlknchn sin 6, = kj (&D cos 6,2 + j, (kr) cos So 
Brun 


Ma (RE) SiN So? /L FLIA cos Se - Notkr) сіп 5212 lom ἘΜ 


For } = 1 

м | = frL cos s, (XRD -2 jatko) - h. si^ &, (Ge) 7 2-08 P)) ]- 
ñ 
risem 


+ CoS S, J, (RO -SIN пин / г (соз (0) - ып 5, Ты БЕ ЕМ 


Rearranging and taking the inverse tangent we get the following express- 


ions for the phase shift: 






сенат! Ва (emt Rick] rT — do кәй - 
[Ze rocked tk ChAT - NORD 


` 


S ES -ва \ (ту - + NC r +y (Rr 


$ = ой (11) 
E n (kO - Вп.) - а n, (Dr +n en 


The total cross section as a function of energy is obtained from 


the phase shift using the relation [ 8] 


2 
cle) = & Ж (191) 51η Sg (12) 
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For identical particle scattering, equation (|?) becomes, for the 


singlet state: 


ote) = $. $Ñ Cans) sin’ sg 


even 


and for the triplet state: 
2 
9, (Е) ъ £ qq em SIM 5; 
odd 


The computation of cross sections was made for both the singlet 
d 20 state and the triplet À = 1 state. The Ñ z I cross sections 

were observed to be at least 3 orders of magnitude smaller in the 
energy range under consideration and were disregarded in subsequent 
computations. The computer program in appendix Ill! however, retains 
the capability of Ñ = 1 calculations. 

The computer program in appendix lll! was used in cross section 
calculations for all the potential parameters in section 3. First 
the program computed the potential parameters at the base scattering 
length and effective range and at percentage deviations from this base 
read into it in its input section. The total scattering cross section 
was then computed as a function of center of mass energy. The calcul- 
ated cross section for each set of potential parameters was then divided 
by the cross section for the base parameters to normalize the cross 
section and present it as a deviation from the base. 

The resulting data was plotted to show the percentage difference 
in total cross section from the calculated total cross section for n-p 
scattering wi th κ... - 23.677 fm and ros = 2.46 fm as a function of 
energy for percentage deviations of a and P The graphs are drawn 
over two distinct energy ranges, 20 - 100 kev and .1-2.6 MeV, to 


better portray results and separate graphs are provided for each 
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potential shape. 


In Moravesik's paper, the function Q(T) is defined by: 


те -ἱ 1 
аст) (44 =+ (0-e = 0/28) 


where Es MT for T s= laboratory kinetic energy and m = neutron mass. 
Moravcsik A determined that in a 10% experiment, a linear approxi- 
mation for Q(T) would be sufficient. 

The value of this straight line at € = 0 would give 1/.2 directly 
while the slope of the straight line would determine Pola and the sign 
of a. 

Moravcsik concludes that a determination of total cross section 
to 10% accuracy at ten suitably distributed energy points could det- 
ermine the scattering length to 3% and the effective range to 50 - 70%. 

Tables IV, V of this report show that a 10% determination of 
total cross section in the vicinity of the n-p scattering parameters 
will allow a 6% determination of ag at 20 kev and a 16% determination 
at 400 kev while at 2.6 MeV a 16% error in scattering length has only 
a 2% effect on total cross section. The same tables show that a 10% 
determination of total cross section will allow a 40% determination 
in effective range at 2.6MeV while a 100% error in effective range 
has approximately a 2% effect at 20 kev. 

An alternate method of determining the effective range parameters 
is therefore offered. The magnitude of the total cross section at 
20 kev could be compared to supplemented data of this report to arrive 
at a value for a, within 64. Graph 17 is shown as an illustration of 
possible supplemental data. The magnitude of the total cross section 


at 2MeV or greater could, in the same manner, be used to arrive at a 
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EFFECT OF A 10% ERROR IN TOTAL CROSS SECTION UPON THE ERROR 


ENERGY (kev) 
20 
60 
100 


200 


% Error in as for 


SQUARE WELL 
8-16 
8-16 
8-16 
8-16 


TABLE IV 


GAUSSIAN 
6 
8.5 

1! 


16.5 


32 


EXPONENTIAL 
6 
8.5 

12 


16.5 


IN a, 
YUKAWA 
6.5 
9.5 
16 
18 








TABLE V 


EFFECT OF A 10% ERROR IN TOTAL CROSS SECTIONS UPON THE ERROR IN e 


Феггог іп ή, for 
᾽ς 
ENERGY (MeV) GAUSS IAN EXPONENTIAL YUKAWA 
1.6 60 58 62 
1.8 56 52 59 
2.0 53 50 57 
2,2 50 47 55 
2.4 45 45 52 
2.6 40 42 50 
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value for r within 50%. Because of the negligible effect of both 


ο,5 


large errors inr at the low energy cross section and large errors 


ο,5 
іп ас at the high energy cross section, they may be determined indepen- 
dently, each having little effect upon the determination of the other. 
Either gaussian, exponential, or yukawa data could be used because of 
only about a 1% difference in a, value among them for any given cross 
section. Thus the scattering length and effective range could be det- 
ermined with a minimum of two points, one at either end of the energy 
range. For greater accuracy, the computer program in the appendix 
could easily be adapted to a searching interpolation. An advantage 
could be gained if the necessity for only two values of cross section 
would allow a more precise determination of those two values. It can 
be seen from accompanying graphs that the energy must be accurate to 
10% to achieve meaningful results in this way. If this energy accuracy 
is not feasible, then the method outlined by Moravcsik for scattering 
length and effective range determination must be used. In any event, 


Moravcsik's development of expected accuracies is confirmed by this 


phase shift analysis. 
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5. Shape independent approximation 


The total cross section for neutron-neutron scattering was cal- 
culated by the shape independent approximation of effective range 


scattering theory using "os and a*_. This cross section was compared 


5 
with that from the potential scattering calculations in section 4. 


7 
Segre gives for the total cross section 


0)E) = 4 | пот TT 
e | +COT ® So gt + Le 


where a(k) is the scattering length as a function of wave number. An 


approximation for a(k) is given for k = 0 as: 


Lc uU uu 2 


) - - zx 
оо) eg seme. ALO» е 


which, when substituted for a(k) in the cross section formula yields: 


4T 


eie» = rae „у ЖК 
2. 


Formula (la) was used in a calculation of total cross section in 
the energy range .l - 2.6MeV and compared with the cross sections 
obtained from potential scattering for the same parameters. The results 
are shown in graphs 32 - 33. 

The results of the effective range theory calculations compare 
favorably with potential scattering calculations. A 2% difference is 
the greatest deviation noted in the comparison. This confirms the 
commonly held belief that the shape independent approximation of 


effective range theory is quite accurate in this energy range. 
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6. Summary 


In section 4 it was determined that a 10% scattering experiment 
could define a scattering length to approximately 6% and an effective 
range to approximately 40% accuracy. This accuracy might be achieved 
with oniy two values of cross section; that at 20 kev and that at 2.6 
MeV, if the energy could be determined to 10% in both regions. Certainly 
it would be achieved by a least squares fit at 10 - 20 appropriately 
distributed energies as proposed by Moravcsik. 

In the region of the accepted n-p scattering length and effective 
range, it was found that a 2% change in scattering length resulted in 
a 3.1% change in cross section at 20 kev center of mass energy while 
only causing a .1% change in total cross section at 2.6MeV. A 20% 
change in effective range was seen to cause a 5% change in cross section 
at 2.6MeV while having almost no effect at 20 kev. Thus the conclusion 
was reached that effective range and scattering length could be determined 
separately and independently by the cross sections at the high and low 
ends of the proposed energy range. 

In section 3 it was found that a 1% change in potential depth 
or range resulted in about a 20% change in scattering length and about 
a 1% change in effective range for Exponential and Gaussian weils. 

For the Yukawa well, a 1% change in potential depth and range 
resulted in about a 12% change in scattering length and a 1% change 
in effective range. 

It was conversely determined that a 6% determination of scattering 
length can determine the Gaussian and Exponential potential depths to 
less than - .2% and the potential range to less than - .3% for the 


n-p scattering length and effective range vicinity. The Yukawa well 
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depth would be determined to less than -.5% and the potential range 
to -.2%. 

It therefore seems quite likely that a 10% scattering experiment 
would determine the scattering length and effective range accurately 
enough to furnish an acceptable set of potential parameters. These 
potential parameters appear to be adequate to investigate charge - 
independent and charge-symmetry effects. 

Effective range theory was quite adequate to determine cross 
sections in the examined energy range. A 1% agreement was found with 
potential scattering. There should be little question of ts propriety 


here. 
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APPENDIX | 


COMPUTER PROGRAM DESCRIPTION 


The program is written in FORTRAN J? language and, as listed in 
Appendix Ill, is arranged to operate with the Control Data Corporation- 
1604 as a central processor and the | B M = 1401 а5 peripheral equipment. 
The library tape functions used are the standard FORTRAN 1! functions 
with the exception of a subroutine called DRAW. This subroutine is a 
part of the library tape at the Computer Center of the United States 
Naval Postgraduate School and is identified by its large computer 
operators co-op 10: J7-NPS-DRAW, The subroutine was written by J.R. 

Ward of the U.S. Naval Postgraduate School. Where this subroutine is 
not available, the program can be run by omitting subroutines SKETCH 
and NOW and their call cards. 

The program consists of a main routine and !3 subroutines. The 
main routine calls the subroutines in proper sequence and also provides 
for a sequential variation of the parameters г,» ала Vo: The subroutines 
and the function each performs are listed as follows: 

Subroutine | МРИТ2 Reads in data cards 

Subroutine READY2 Calculates and tabulates energies, wave numbers, 
scattering lengths, effective ranges, potential 
parameters, and calls bessel functions. 

Subroutine FIGURE Prepares for and calls integration routine, forms 


the logrithmic derivative of the integrated function. 


Subroutine RKUTTA Integrates equation provided by DERIV for f s O. 
Subroutine RKUTTAI Integrates equation provided by DERIV] for f = 1. 
Subroutine DERIV Provides Qs O equation for integration for a 


particular potential. 
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Subroutine 


Subroutine 


Subroutine 


Subroutine 


Subroutine 


Subroutine 


Subroutine 


DERIVI 


FINISH 


OUTPUT 


BESSEL 


SKETCH 


APPEND 


NOW 


Provides 1- | equation for integration for a 
particular potential. 

Computes phase shifts, cross sections, and 
normalized cross sections. 

Prints out resulting ro/a ratio, Vo, ra, energy, 
phase shifts, and cross sections. 

Calculates bessel functions 

Prepares for and calls DRAW to plot all! cross 
sections on one graph for comparison. 

Stores base cross section for normalization 

in FINISH, 

Prepares for and calls DRAW to plot all norma- 


lized cross sections on one graph for comparison. 
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EINT 


E I NC 


ROS 
AS 


CHANGE 


MANY 


REND 


NUMBER 


MODE 


MOMENTA 


FK 
FK2 
ЕМАКУ 


КА 


APPENDIX I! 


PROGRAM DICTIONARY 


Initial value of center of mass energy in MeV. 
Incremental value of center of mass energy in MeV. 
Maximum number of center of mass energy values. 
Base effective range in cm. 

Base scattering length in cm. 

Percentage incremental change in base effective range 
or base scattering length in decimal form. 

Number of + or = incremental charges to be made. 
Radius at which integrated solution is matched to 
exact solution for phase shift determination in cm. 
Number of integration steps. 

Print out control for program monitoring 

а 0 {ог result print out only 

= | for minor progress report 

- 2 for print out of all calculations 

Energy subscript for wave function print out. 

20 Xs 0 calculations only 

zl = 1 calculations only 


2 both E O and 21 calculations. 


Wave number, k 

Wave number squared, к^. 

MANY in floating point. 

Ratio of effective range to scattering length, rj/a. 


Well depth parameter, 5 
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РОТ 

КМ 

АКСВ 
VALUE 
VAL JO 
VAL JI 
VAL J2 
VALNO 
VALNI 
VALN2 

A 

V 
FNUMBER 
STEP 

R 

YVARS (1) 
YVARS (2) 
XVAR 

Y 


YDOT 
COMPARO 
COMPAR | 


XC 


Intrinsic range, b 

Well depth, Vo 

Well range, rn 

Bessel function argument, kr 

Output array of bessel functions 

Spherical Q = 0 bessel function, j 

Spherical = 1 bessel function, j| 

Spherical Qs 2 bessel function, jo 

Spherical {= 0 neumman function, по 

Spherical > 1 neumann function, п, 

Spherical ἢ » 2 neumann function, n2 

A particular potential range from array, RN 

A particular potential depth from array, POT 
Number of integration steps in floating point 
integration step size 

Array of radius values at each integration step 
Wave function 

Slope of wave functions 

Radius variable of integration 

Wave function at end value of integration step for f= 0. 
Slope of wave function at end value of integration 
step for (= o. 

Logrithmic derivative of wave function at end of 
integration range for ( s o. 

Logrithmic derivative of wave function at end of 
integration range for P= 1. 


Runge kutta weighted radius step 
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YC 

DY 

YL 
YLDOT 
DELTAO 
DELTAI 
XSECO 
XSECI 
TOTAL 
BETTER 
FNORM 


NUMBERR 


Dependent variable for XC 
Derivative of dependent variable for XC 
Same as Y for 1- 1. 
Same as YDOT for } = ja 
Y = 0 phase shift. 
$ = 1 phase shift. 
À = O cross section 
| = | cross section 
Sum of singlet and triplet cross sections 
Normalized {= O cross sections 
Base cross section to be used for normalization 
Number of simultaneous first order differential equa= 


tions to be integrated. 
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APPENDIX III. 
PROGRAM LISTING 














MAIN PROGRAM 

DIMENSION E(100)sFK(100) »FK2(100) »RN(100) »VALJO(100) sVALJ1( 100) 
1»VALJ2(100)» VALNO(100), VALN1(100)5»VALN2(100)» R(1000), Y(1000), 
2YDOT(1000), COMPARO(100), COMPAR1(100), YL(1000),YLDOT(1000),POT(1 
300)» DELTAO(100) )DELTA1(100),XSECO(100) »XSEC1(100)  TOTAL(100) „КА (10 
40) sFNORM( 100) , BETTER (10,100) 

COMMON E»SEINT»EINC, FK, FK2» RINT, RINC», RN, A» 11» LL» NN» MODE, 
1REND>VALJO»VALJ1)VALJ2>VALNO>VALN1>VALN2>JJ3M, NUMBER» R>I»L>V>K> 
2N>Y» YDOT» YL» YLDOT» Ill» COMPARO» COMPAR1>POT>VINT>VINC>+MOMENTA» 
ЗОЕТАО „ОЕ| ТА1 , Х5ЗЕСО , ХЗЕС1 „ ТОТА. „ВА „ВО5 »А5 CHANGE >) MANY S.B, FNORM,BE 
4TTER 

CALL INPUT2 

CALL READY2 · 

A=RN(MANY+1) 

V=POT (MANY+1)/.41499E-24 

CALL FIGURE 

CALL FINISH 

CALL APPEND 

DO 1 ритми 

A=RN(L) 

= POT(L)/+41499E-24 

CALL FIGURE | | 

CALL ‘FINISH 

CALL OUTPUT2 

CALL SKETCH(E»XSECO»MODsNUMPTSo IIs Lo LL! 

CONTINUE 

CALL NOW(E>BETTER»II>LoLL) 

END 


ον 
RÀ 





SUBROUTINE INPUT2 

DIMENSION E(100) »FK( 100) »9FK2(100) »RN( 100) sVALJO(100) »VALJ1(100) 

1 „МА .)2(100)» МА!МО(100)» МА М1 (100) „МАСАЖ (100), К( 1000)» У( 1000)» 
2YDOT(1000), СОМРАКО( 100)» СОМРАКТ (100)» У| (1000) „УГРОТ (1000) »РОТ(1 
300) sDELTAO(100) »DELTA1(100) »XSECO( 100) »sXSEC1 (100) sTOTAL (100) »КА (10 
40) »FNORM(100) sBETTER (109100) | 
COMMON ESEINTSEINC» FK» FK2» RINT» RINC, RN, A, Il, LL, NN, MODE, 
1REND»VALJOS,SVALJ1»VALJ2s»VALNOsSVALN1»VALN25» JJ, M, NUMBER» R»Is5»L5,V5K, 
2N»Y» YDOT, YL» YLDOT» III» COMPARO» COMPARI»POT»,VINT»VINC»MOMENTA » 
3DELTAO»DELTA1»XSECO,XSEC1»TOTAL,RA,ROS,AS,CHANGE,MANY»S,B, FNORM,BE 
4 ТТЕК 

READ 1»EINT>EINC>1I7 

ЕОКМАТ ( 2Е20.4» 110) 

‚READ 2,ROS»AS»CHANGE »MANY 

:FORMAT(2E20.6,F10.4, 110) ( 

READ 39REND» NUMBER 

FORMAT(E20+.49+110) 

READ 4»MODE. 

.ҒОВМАТ( 110) ` ! 

BREAD 5», Ill I ; 

FORMAT ( 110) ` 

READ 6, МОМЕМТА г 
FORMAT (110) 

IF (MODE ~ 1) 858,7 

РКІМТ 1ЭЕІКТ»ЕІКС»ІІ 

PRINT 2,ROS,AS»CHÄNGE ¿MANY 

‘PRINT 39RENO NUMBER; 

PRINT 4 sMODE 
PRINT 5,111 
‚PRINT 6»MOMENTA 
RETURN 

END 


“м. 
4 “~ 
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SUBROUTINE FIGURE 

DIMENSION E( 100) sFK(100) sFK2(1900) sRN( 100) ΝΑΙ ώΟ(100) »VALJ1(100) 
1sVALJ2(100) 5 VALNO(100)s VALN1( 100) sVALN2(100)s R(1000)s Y(1000)» 
2Ү00Т7(1000)» СОМРАКО(100)» СОМРАК1(100)» У1 (1000) „УГООТ 1000) „РОТ: 1 
300) sDELTA0(100) sDELTA1 (100) »Х5ЕСО (100) »Х5ЕСТ (100) +» ТОТАС (100) »ВА (10 
40) „ЕМОКМ( 100) sBETTER(109100) 

DIMENSION YVARS(30) 

COMMON EsSEINTsEINC» FKs FK2» RINT» RINCs» RN, As» II» LL» NN, MODE» 
I1REND»VALJO,VALJ1»VALJ2»VALNO »VALNI1»VALN2 »JJ,M, NUMBER, R»I;L,V»K; 
2N>»Y» YDOT»> YL» YLDOT» III» COMPARO> COMPAR1>POT>VINT>VINC9MOMENTA > 
ЗОЕСТАО „ОЕГТА1 ,Х5ЕСО , ХЗЕС1 ,» ТОТА: „КА „КО5 »А5, СНАЛОЕ » шы В,ЕМОКМ › ВЕ 
4TTER 
FNUMBER = NUMBER 

STEP = REND / FNUMBER 


R(1) = O, 
NN = NUMBER + 1 n 
: DO 299 М<-2, „ММ l 
299 R(N) = R(N-1) + STEP 
DO 292 К«2,11 
"YVARS(1) - 0. 
YVARS(2) = le 


IF (MOMENTA - 1) 283532845283 
283 CALL RKUTTA(2»XVAR>»YVARS»STEP) 
IF(K-III) 29392979293 | 
297 PRINT 296 
296 FORMAT (13X5,7H ENERGY 513X»7H RADIUS»6X»14H WAVE FUNCT ION) 
PRINT 295,E(K) 
295 FORMAT(F20.2) 
PRINT 2945 (R(N)s YIN)» N&loNN) 
294 FORMAT (E40.4>E20+4) 
293 COMPARO(K)=YDOT(NN=1)/Y(NN=1) 
284 IF (MOMENTA - 1) 29252915291 
291 YVARS(1) - 0. | 
YVARS(2) = le 
CALL.RKUTTA1(2>»XVAR»YVARS»STEP) 
IF(K-III) 28592899285 
289 PRINT 288 
288 FORMAT (13Xs7H ENERGY 513X%7H RADIUS»6Xs18H WAVE FUNCTION L= 1) 
PRINT 287,E(K) 
287 FORMAT(F20+2) 
PRINT 286» (R(N)+YL(N) +N=1,NN) 
286 FORMAT(E40.4>E20.4) 
285 .COMPAR1(K)-YLDOT(NN- 1)/YL(NN-1). 
292 CONTINUE 
RETURN 
END 


کس سک د ا --+< - - - — —— о‏ — 
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SUBROUTINE RKUTTA(NUMBERR sXVARsYVARS»sSTEP) 

DIMENSION E(100) „ЕК: 100) „ЕК (100) »RN( 100) sVALIJO(100) »VALJ1(100) 
1ЭУА(/2(100)» УА(М0(100)» VALN1(100)»VALN2(100)» R(1000)5, Y(1000)» 
2YDOT(1000)» COMPARO(100)» COMPAR1(100)» YL(1000)»YLDOT(1000)»POT(1 
300) sDELTAO( 100) sDELTA1(100) »XSECO(100) »XSEC1(100) sTOTAL (100) sRA(10 
40) sFNORM( 100) »BETTER(109100) 

DIMENSION YVARS(30) sAK (4930) sDY(30) sYC(30) 9C (4) 

COMMON EsEINTsEINCs FKs FK2s RINTs RINC»; RN, A» II» LU» NN» MODE» 
LREND »VALJO » VALJ1 » VALJ2 » VALNO sVALN1 sVALN2 9JJ9My NUMBER» RoIoloVoKsy 
2М,Ү» ҮбОТ» ҮС» ҮШбОТ» ІІІ» СОМРАВО,») COMPAR1>POT>VINT>VINC>MOMENTA» 
3DELTAO sDELTA1 »XSECOsXSEC1sTOTAL »RA „КО5 »А5 „СНАМОЕ, МАМУ » 5 „В, Е МОКМ „ВЕ 
4TTER i 

C(1)=0.0 

C(2)20,5 

C(3)=0.5 

C(4)=1.0 

DO 298 ΝΞΙΝΝ 

XVAR = RIN) | 

ШІ! 1-1,4 0 

XC=XVAR + C(I)*STEP 

DO 2 J=1»NUMBERR Ἢ ` 

ҮС(Ј) =YVARS(J) + C(I)XAK(I-19) 

BEALL DERIV(XCsYCsDY) 

DO 1 J=l»NUMBERR 
ИКТ, 1) = STEP * DY(J) 

DO 3 J = 1» NUMBERR. : 

YVARS(J) zYVARSUJ) * (AK(1»J) £269 *XAK (232) 29 *AK(3»3 J) +ҒАК( 497) ) /6е 


PYDOT(N) = YVARS(2) ~ 
-Y(N) = YVARS(1) 

RETURN 

END 








SUBROUTINE RKUTTA1 (NUMBERR sXVAR»sYVARS » 5 ТЕР ) 

DIMENSION Е(100) „ЕК( 100) „ЕК? (100) „КМ (100) „МАГ 10 (100) „МА: .)1 (100) 

1 „МА: .121100) » МА МО( 100)» МА М1СТО0) УЧАС М2 (100) 5 σος ο. 
2YDOT(1000)» СОМРАКО( 100)» СОМРАК1 (100), У1 (1000) „УГРОТ (1000) „РОТ (1 
300) „ВЕ! ТАО (100) „ВЕ! ТА: (100) „ ХЗЕСО (100) „Х5ЕС1 (100) „»ТОТА| (100) „КА (1С 
40) „ ЕМОКМ (100) „ВЕТТЕВ (10»100) 

'DIMENSION YVARS(30)sAK (4930) „ОУ(30),УС(30) „С (4) i 
COMMON ESEINT»EINC» FK» FK25, RINT» RINC», RNs A» II» ШІ» NN» MODE; 
IRENDS,VALJO,VALJ1.5VALJ25VALNOSVALN15VALN25 JJSM, NUMBER, R5»Is»L5,Vs5K; 
2N>»Y» YDOT» YL» YLDOT>» III» COMPARO, COMPARI1SPOTS,VINTsSVINC»MOMENTA 
3DELTAO,DELTA1»XSECO,XSECI1» TOTAL ,RA,ROS S5 AS, CHANGE, ΜΑΝΥ 5280 ЕМОКМ»ВЕ 

4TTER 


C(1)=0.0 
C(2)20,5 
C(3)=005 
C(4)=1.0 


DO 290 N=1»NN 
XVAR = К(М) 
EDO l I=1,4 | | 
'XC=XVAR + C(1)*STEP | ; 
¿DO 2 J=1»NUMBERR .. ` | 
2: YC(J] =YVARS(J) + C(1)*AK( 1-19.) 
¡CALL DERIV1I(XCsYCsDY) 
г ВО 1 J=1»NUMBERR 
l'AK(I»sJ) с 5ТЕР * DY (J) 
DO 3 J = 1yNUMBERR 
3:YVARS(J) sYVARS(J) + (AK( 19d) +2e*AK( 295) 4+20*AK( 395) +AK(49J5))/60 


СҮМ) = YVARS(1) 
290; YLDOT(N) = YVARS(2) 
RETURN 
END 
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SUBROUTINE ЕТ НТ 

DIMENSION Е(100),ЕК( 100) „ЕК? (100) »RN( 100) я МАЦ ЈО (100) „МА 21 1100) 

1 » МА: .)2(100)» МА МО (100) » МА| М1 (100) „МАГМ2 100)» К( 1000)» Y(1000)» 
2YDOT(1000)» СОМРАКО (100), СОМРАК1 (100) » МУ. (1000) „УГООТ (1000) „РОТ (1 
300) sDELTAO( 100) sDELTA1 (100) »XSECO0(100) »XSEC1(100) sTOTAL (100) sRA(10 
40) „ЕМОКМ (100) ЭВЕТТЕК(10»100) 

COMMON E,.EINT»EINC, FK, FK2» RINT»; RINC, RN, A» II» LL» NN», MODE, 
]J)REND»VALJO»VALJ1»VALJ2»VALNO»VALNLI5»VALN25»JJs»M» NUMBER» R»5»Is5»L5V5K» 
2N>Y» YDOT» YL» YLDOT» III» COMPARO» COMPAR1>POT>VINT>VINC>MOMENTA » 
3DELTAO»DELTAlLl»XSECO»sXSEC1»TOTAL»RA,ROS,AS»CHANGE ,»MANY,S,B» FNORM ΒΕ 
ENTER 

IF(MOMENTA - 1) 39953979399 


ESO DO 398 Із2,11 


298 
97 
396 


095 
290 
569 
394 


900 


DELTAO( 1) =ATANF ( (REND*® (COMPARO( 1) ¥VALJO( 1) +FK(1)*VALJ1(1))= 
1VALJO(1)) /(REND*(COMPARO(1)X*VALNO( I)+FK(I)VALN1(1))- 
2VALNO(1))) 

XSECO(1)=4.*3.1416/FK2(I)*SINF(DELTAO(1))**2 

ТЕ ( МОМЕМТА- 1) 394 396396 

DO 395 12, 11 

ОЕ: ТА! (Г) <АТАМЕ ( ( КЕМОХ (ЕК (Т) /З. “МУАС.)О(Т)-2./3. ЖЕК(Т) Х 
1МАЦ.)2 (ГТ) -СОМРАВ 1 (Г) МАС.) 1 (Г) ) + МАС.1 (Г) ) / (ВЕМОЖ (ЕК (Т) /З. 


2УА|МО(Г) -2. /3.ЖЕК (Г) КУАМ2 (Г) -СОМРАВТ (Т) “МА М1 (Т)) + МА М1 (Т))) 


XSEC1(1)=12.*3.1416/FK2(I)*SINF(DELTA1(1))**2 

DO 389 I - 2,11 - 

TOTAL(I) = «25 * XSECO(I) + «75 # Х5ЕС1(Т) : 
DO 900 I-22,II 1 

BETTER(L»I) s XSECOLI)/FNORM(I) 

RETURN | 

END | 
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SUBROUTINE OUTPUT2 
DIMENSION E(100)»FK(100)»FK2(100)»RN(100)»VALJO(100) »VALJ1(100) 
11, VALJ2(100)» VALNO(100)» VALN1(100),VALM2(100)» R(1000)» Y(1000)» 
2ҮрОТ (1000), СОМРАКО (100), СОМРАКТ( 100)» У1 (1000) „УГРрОТ (1000) „РОТ (1 
300) „ВЕ | ТАО (100) „ОЕ ТА (100) +» ХЅЕСО (100) » ХЅЕСІ1 (100) + ТОТАС (100) » ВА (10 
40) „ЕМОКМ( 100)» ВЕТТЕК (10»100) | 
СОММОМ ЕЭЕІКТ»ЕІКС» ҒК» ЕК2› КЇМТ› КЇМС, КМ, А› 11» LL» NN» MODE, 
)REND»VALJOsVALJ1sVALJ2sVALNO»VALN15»sVALN2»JJ»M» NUMBER» R»I»L»V»K» 
2N,Y» YDOT, YL» YLDOT» III» COMPARO, COMPAR1,POTS,VINT5VINC$»MOMENTA » 
3DELTAO sDELTA1 »XSECOsXSEC1s TOTAL sRA »sROSsASs CHANGE sMANY »S5BsF NORM, BE 
4TTER 
PRINT 600,RA(L) 
600 FORMAT( 34H EFFECTIVE RANGE/SCATTERING LENGTHs/3E34e10) 
PRINT 599,A,POT(L) 
599 FORMAT (13Xs7H RADIUSs4X916H POTENTIAL DEPTH s/ sE20¢4sE20 04) 
PRINT 598 
598 FORMAT( 7H ENERGYs12H PHASE ANGLE s15H CROSS SECTION» 12H: PHASE ANGLE 
l»15H CROSS SECTIONs2X96H TOTAL s/s11Xs4H L=097X94H LEO, Ον. Lal, 
28X 94H L=1s8Xs6H CROSS»/»63Xs8H SECTION) 
PRINT 597, (ELI) SDELTAO LI) »XSECOCI 3DELTAL CI) 9XSECL CI Io TOTAL CI] o Ds 2 
0511) 
597 FORMAT(F7.2,E12.4,E15.4,E12.4,E15.45 Elle 4) 
PRINT 596 
596: FORMAT(7H ENERGY s5Xs25H NORMALIZED CROSS SECTION} 
PRINT 595,(Е(І)»ЭВЕТТЕК((»1>1І2,11) ~. 
595 ЕОВМАТ(Ғ7.2>Е20.8). | 
RETURN 
END 





s OO 


299 
298 


300 


10 


SUBROUTINE BESSEL (Х,М,МАГОЕ,ТМОЕХ,ММ) 
DIMENSION VALUE (40) 5B4(100)sF (100) 
IF (X) 1091099 


ІК-Х/2.0 

K=N+20+IK 

IF(K-100) 30092995299 

PRINT 298 

FORMAT ( 66H ARGUMENT PLUS ORDER TOO LARGEe INCREASE DIMENSIONS AND 


1RUN AGAIN.) 


RETURN 

CONTINUE 

2-1.0/Х 

РІ-3.141592653 
GOTO(405405605»605,5805905,51005110) » INDEX 
PRINT 8,X 


80FORMAT (1X»63HSUBROUTINE BESSEL DOES NOT WORK FOR X0 OR LESS THAN 


40 
Ll 
E13 
ἘΠ 
14 


43 
45 


42 
ma 


20 
24 


15 


1. 0» HERE X= £15.48) 


RETURN \ 
BJ(K+1)=0+0 : 

BJ(K)=10.0E-30 К 

L=K s 
FL=L | 
BJ(L-1)22.0*FL*ZX*BJ(L)-BJ(L-*1) 

Me(t—-1) 12912913 N 

П-(-і = 
8o TO 11 

SUM=0.0 

IK=K 

DO 14 I=2»IK»2 

SUM=SUM+BJ(1) 

Cz1.0/(BJ(0)-424 OX* SUM) 


GO TO(425423)5» INDEX , 


DO 45 I=09»K 

BJ(1)=C*BJ(1) 

GO TO 54 

DO 41 І=0,М 

VALUE(1)=C*BJ(1) 

RETURN 

GO TO 40 

EC=.5772156649 

SUM=0.0 

IK=K/2 

рО 15 Ізізік 

Е1< 1 | 
50М-50М%((-1.0)%%( 1-1) %В./(2%1))/ҒІ 
Е(0)-(2.О/РІ) %(ВӘУ(03%((ОСЕ(Х/2.0) «ЕС)%2.0%51)М) 


F(1)=(BJ(1)*XF(0)-2.0/(PI*X))/BJ(O0) 


IN=N-1 
DO 16 Is1»IN 


EI-I 





16 F(I*1)224,0X*FI*ZX*F(I)-F(CI-1) 
DO 51 I=0»N 
51 VALUE(I)=F(I) 
RETURN 
60'BJ(K+1)=00 
BJ(K)=10.E-30 
“(аҚ | 
19 Ғісі 4 : 
BJ(L-1)22,. *FL*XZ*BJ(L)-BJ(L*1) 
IFtE=131 75172123 
18 L=L-1 
EGO ΤΟ 19 
17 SUM=0.0 
сро 21 151. 
21 SUM=SUM+BJ(I1) 
„С (ЕХРЕ (Х))/ (В.) (0) +2 Ж5ЦМ) 
GO TO(15»1561562) INDEX 
ІШ РВІМТ 1000 
=- RETURN 
1000 FORMAT(40H CARD MISSINGe IGNORE RESULTS. 
62:00 64 1:0,1 
64 BJ(1)=C*BJ(1) 
moo TO 71 
61 DO 65 [=0,М 
65 'VALUE(I)=C%*BJ(I) 
RETURN 
MO. GO TO 60 | 
71 1Е(5.0-Х) 120,121,121 
121 M=10 Dio 
60 TO 122 
120 M=5 
122 F (0)=GAUSS(0e091e09MsX) 
F(1)=(1./X-F(0)*BJ(1))/BJ(0) 
ІМ-М-1 | Е 
DO 22 I=lsIN 
ΕΙΞΙ : 
22 F (I*1)5224,0*FI*Z*F (ІУ)ЖҒ (1-1) 
ро 72 ΙΞ0.Ν 
V2 VALUE(I)=F (I): 
RETURN 
80 Е(К+1) +0.0 
F (K)=10.0E-30 
„ЕК 
25 FL=L 
F (L-1)=(2.0*FL+1.0)*2*F (L)-F (L+1) 
IF (L-1) 23,23,24 
24 L=L-1 
КОО το 25 
23 C=SINF(X)/(X%F (0)) 
ТЕ (ММ) 200,200»201 
201  SQ=SQRTF(2.*#X/PI) 
DO 202 ΙΞ0.Ν 
202 :VALUE(I)=C#SQ#F (I) 
RETURN 


m 








200 DO 26 І-0»М 
26 VALUE(I)=C*F (1) 
RETURN | 
90 F (0)5-COSF'X) 7X 
F (1)2-SINF(X)/X-COSF(X)/X**2 
IN=N-1 
DO 27 Іс191М 
ВЕТ 
27 F(I+1)=(2.*FI+1.0)%*ZXF(1)-F(I-1) 
ТЕ (ММ) 20332035204 
204 SQ=SQRTF(2.*X/PI) 
DO 205 I=09+N 
205 VALUE(I)=F (1)%50%(-1.)%%(-1-1) 
RETURN 
E203.DO 91 1=0›М 
91 VALUE(I)=F (І) 
‘RETURN 
ШОО E (K+1)=0.0 
F (K)=10.0E-30 
L=K 
30 FL=L 
I (L-1)=(2+sO%#FL+1,0)%Z*F (LI+F (L+1) 
IF (L=1) 28928529 
29.L=L-1 
GO TO 30 
28 C=SINH (X)/(X*F (0)) тт. 
IF(MM) 20692063207 | 
207 SQ=SQRTF(2+*X/PI) 
DO 208 ΙΞ0.Ν 
208 VALUE(I)=C*SQXF (І) 
| RETURN 
206 DO 101 I=OyN 
101 VALUE(I)=C*#F (I) 
RETURN 
110 Е (0)--РІЗЕХРЕ(-Х)/(2.0%Х) 
Е (1)+ (Р1# (Х+1.0) ЖЕХРЕ (-Х)) / (2 .„Ожхжж2) 
ІМ-М-1 
DO 31 ]=1,IN 
FIs] | 
31 F (1+1)=F (I-1)-(2.0*FI+1.0)*Z*F (1) 
IF(MM) 20952093210 
210 50-5ОВТҒ(2.%Х/РІ) 
МО 211 Іс0,М 
211 УАГМЕ(І)-Ғ (1)%50%(-1.)%%(-1-1) 
RETURN 
ӘС СОО 111 1=0,м 
111 VALUE(I)=F (1) 
RETURN | 
ЕМО š р 
` FUNCTION GAUSS NUMERICAL INTERGRATION 
; FUNCTION GAUSS REQUIRES FUNCTION GAUSSF(X) NUMERICAL INTERGRATION 
FUNCTION GAUSS(G»B»M»X) 
.DIMENSION R(10)5U(10) 
‚Ut1)=+.0744371695 








0(2)--.0744371695 
U(3)=+02166976971 


0(8)-.2166976971 
U(5)=+.3397047841 
U(6)=-.3397047841 
U(7)=+, 6325316833 
U(8)=-,&325316833 
U(9)=+.4869532643 
1110) +-,4869532643 
R( 1)=+. 1477621124 
R(211147 62120 

‘R(3) =+e1346333597 | 
R(4)=++.1346333597 
В(5) =, 10954318113 
В(6) =, 10954318113 
R(7)=+.07472567458 
R(8)=+.07472567458 
R(9)=+.03333567215 
R(10)244.03333567215 
A=G 
FN=M 
P=(B-A)/FN | 
GAUSS=0.0 
DO 1 J=1»M 
С-0.0 
Y=A+P 
DO 2 I=1»10 
D=(R( 1) *GAUSSF( (Y=A) ¥U(T)+(A+Y) /2009X)) ¥(Y-A) 
2 C=D+C i 
ΑΞΥ 
1 GAUSS=GAUSS +C 
RETURN 
END 
FUNCTION GAUSSF(U»X) 
IF (1e0/U-710.0) 12391245124 
124 GAUSSF=0.0 
RETURN 

123 GAUSSF= (EXPF(-X*COSH(1+0/U-1.0)))/U*%*2 

END | 


qe 








SUBROUTINE 5КЕТСН(Е» ХЅЕСО » МОр » МОЈМРТ5 » ІІ» LL) 
DIMENSION Е (100) » ХЅЕСО (100) » ХЅЕС1 (100) F TOTAL (100) 
DIMENSION X(900)5,Y(900) 5. ITITLE(12) 
LABEL = 4H 

DO 1 I=1>12 

υπ τες) = ӘН 
ITITLE(1)=8HOSGOOD Р 
ITITLE( 2) =8HeBe | BOX 

ITITLE(3)=8HO9 

ITITLE(4) =8HEXPONENT 
EIITLE(5)=8HIAL WELL 

ITITLE(6)=8H RUN 15 

epi TL&(7)=8H GRAPH O 

IN TEE (8)=8HF CROSS 
ШИГЪБЕ (9) ЕЗНЪЕСТТОМ 
ШЕГТІ.Е(10) + 8НУ5 ENERG 
ITITLE(11)=8HY FOR L= 

ITITLE( 12) = 8HO 

NUMPTS = 11-1 

Mat ti-1)2s293 

МОО + 1 

GO TO 6 

σι 59594 

MOD = 2 

GO TO 6 

MOD = 3 

CONTINUE 

CALL DRAW(NUMPTS>»E»XSECO»MOD»0O>»LABEL»ITITLE>.3>»0>0»0>0>0»9»>»15»>1>»L 
LAST) | . 

RETURN 

END 








SUBROUTINE APPEND 

DIMENSION E(100)sFK(100) sFK2(100) »RN(100) sVALJ0(100) sVALJ1(100) 

1 »МА| 22 (100), МА|МО(100), МА М1 (100) sVALN2(100)» R(1000)»s Y(1000)»5 | 
2YDOT(1000) +» COMPARO(100)s COMPAR1(100)»5 YL(1000)sYLDOT(1000) «ΡΟΤ(1 
300) „ВЕ| ТАО (100) „ОЕ| ТА: (100) „Х5ЗЕСО (100) sXSEC1 (100) sTOTAL( 100) » ВА (10 
40) „ЕМОКМ (100) »ВЕТТЕВ (10»100) | | 

COMMON ESEINTSEINC, ЕК» ЕК2, RINT» RINC» RN» As Il» LL» NN, MODE, 
ıREND>VALJO,VALJ1»VALJ2,VALNO»VALNI1»VALN2,JJ,M, NUMBER» Rs; L; V; K; 
2МҮ» ҮрООТ» ҮШ» ҮШрОТ» І1І» СОМРАКО» СОМРАК1»РОТ»ЭУІМТ»ЭМІМС»МОМЕМТА, 
ЗОЕСТАО „ОЕ| ТА! „Х5ЗЕСО , Х5ЕС1 , ТОТА: „ВА , КО5 „А5, СНАМОЕ , МАМУ »5 „В, ЕМОВМ „ВЕ“ 
4TTER 

ро 1 1=2,11 
1 FNORM(I) = XSECO(I) 

RETURN 

END 
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SUOGROUTINE NOW(ESBBITERsSILSISSEL) 
DIMENSION E(100)» ВЕТТЕК( 10»100) »Х (900) sY¥Y(900) sITITLE(12) 
LABEL = 4H 
DO 1 1=1›12 
1 ITITLE(I)=8H 
ITITLE(1)=8HOSGOOD F 
Ш-ІТІЕ(2)-8Н»В.эВОХ 
ITITLE(3)=8H09 
ITITLE(7)=8HNORMALIZ 
BIITLE(8)=38HED CROSS 
ITITLE(9)=8H SECTION 
ITITLE(10)=8H FOR L=0 
NUMPTS = II-1 
DO 7 L=1»LL 
I (L—1) 2253 
2 MOD=1 
GO TO 6 
3 IF(LL—-L) 5,5,4 
& MOD = 2 
GO TO 6 
MOD = 3 
CONTINUE 
DO 20 1-2,11 
ЕО БЕТ(І)-ВЕТТЕРК(С>1) 
CALL DRAW(NUMPTS+E»BET» MOD»0O»LABEL>»ITITLE>.3>0909090»099>15>1>L: 
1AST) 
CONTINUE А 
RETURN - 
END 
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THE FOLLOWING SUBROUTINES ARE DEPENDENT UPON A PARTICULAR POTENTIAL SHAPE. 





99 


170 
197 
196 
195 
E71 
172 
194 
193 
1192 


192 
192 


GAUSSIAN POTENTIAL 


SUBROUTINE READY2 

ОТМЕМЗТОМ Е1100) „ЕКС 100) „ЕК (100) »RN(100) sVALJO(100) ΝΑΙ 1 100) 
1»VALJ2(100),; VALNO(100)» МА|М1 (100) „УАМ2 (100) » К( 1000)» Y(1000)» 
2700Т (1000), СОМРАКО (100), СОМРАКТ (100) » У1 (1000), УГООТ (1000) „РОТ 1 
ВО) „ОЕГТАО100) »ОЕТАТ (100) »Х5ЕСО : 100) »XSEC1 (100) » fOTAL( 100) sRA(10 
40) „ЕМОКМ (100) sBETTER(105,100) 

DIMENSION VALUE (40) 

COMMON ES,.EINTSEINC, FK, FK25» RINT» RINC» RN» As IIs LL», NN, MODE» 
ÄREND»VALJO»sVALJ1»VALJ2»VALNO»VALNI»VALN2»JJ»sM, NUMBER» Ro»oI>L>V>K» 
2N»Y» YDOT» YL»s» YLDOT» IIIs COMPARO, COMPAR1,POT,VINT,VINC,MOMENTA, 
3DELTAO,DELTA1»XSECO,XSEC1»TOTAL»RA,ROSsAS»CHANGE »MANY»S»B» FNORM ,BE 
ТРЕК 


E(1) = EINT 
DO 199 I = 2,1I 
EN 5 Е(1-1) + EINC 


ЕКТ) = 5ЗОКТЕ (Е (Т) /.41499Е-24) 

FK2(I)2 FK(I) * FK(I) 

LL=2*MANY + 1 

JJ=LL 

FMANY = MANY 

GAUSSIAN DISTRIBUTION 

RA(1) = ROS/AS*(1l.- FMANYXCHANGE) 

DO 170 L=1>»LL l 
5:(1.-2.ОЗКВА(|) + .587ХКА (|) #2) /(1.-2.81 “КА(Г) + 1. ВОБУКА (|) ЖК?) 
55=5-1. ` 

В-КА(()) ҰД5/(1.-.641%55-4.568%55%%2-,721%55%Ұ%23-41.172%55%%4-1.015%55% 
1%5) 

POT(L)=S*.229208E-23/B**X2 

RN(L)=SQRTF (B*X*2/2.0604) 

RA(L+1) = RA(L) + CHANGE *ROS/AS 

IF(MODE - 1) 19051925197 

PRINT 196 
FORMAT (13X,7H ЕМЕКСҮ»5Х912Н МАУЕ NUMBER,3X,20H WAVE NUMBER SQUARE 
1D) 
PRINT 195 ,(E(I)» FK(IOSFK2(I)5I225II) 
FORMAT (3E20.4) 
PRINT 171 
FORMAT(34H EFFECTIVE RANGE/SCATTERING LENGTH) 
PRINT 172»(RA(L) L=]l»LL) 

FORMAT(E34.10) 

PRINT 194 

FORMAT (13X»7H RADIUS) 
PRINT 193» (RN(L)» L=leLL) 

FORMAT (Е20.4) 

PRINT 1192 


FORMAT (16H POTENTIAL DEPTH) 
PRINT 2192» (POT(J)»J=19JJ) 
FORMAT (F16e2) 


PRINT 191 
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189 
188 
187 


186 
185 
184 


FORMAT( 34H RADIUS AND ENERGY VALUES COMPUTED) 
Ρο Вт ите 2 Г] 

АКОВ = ЕК(І) * REND 

CALL BESSEL( ARGBs2sVALUE 95s 0) 
VALJO(I) = VALUE(O) 

VALJ1(—(1) w VALUË (1) 

VALJ2(1) = VALUE(2) 

CALL BESSEL (ARGBs2sVALUE» 690) 
VALNOLI) = VALUE(O) 

VALNI( I) = VALUE(1) 

VALN2(1) = VALUE(2) 


IF(MODE — 1) 18791879189 
PRINT 188,1 

FORMAT ( 1415) 

CONTINUE 

IF(MODE - 1) 184,186,186 

PRINT 185 

FORMAT( 26H BESSEL FUNCTIONS COMPUTED) 
RETURN 

END 
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πο 


ET 7 
196 


Ж 
ТА! 
172 
194 


93 


192 


E92 
E92 


EXPONENTIAL POTENTIAL 


SUBROUTINE READY2 

DIMENSION E(100)»FK(100)»FK2(100)»RN(100) »VALJO(100) »VÁLJ1(100) 
Me ALEJ2(100), VALNO(100)> VALN1(100)»VALN2(100)9 R(1000)» Y(1000)% 
2YDOT(1000)» COMPARO(100)» COMPAR1(100)» YL(1000),»YLDOT(1000)»POT(1 
ОТЭОЕЦГТАО(100)ЭбЕСТА1(100)»Х5ЕС0(100)Х5ЕС1(100)»ТОТАС(100)»КА(10 
40) „ЕМОКМ (100) sBETTER(109100) 

DIMENSION VALUE (40) ; 

COMMON ЕЭЕІМТ»ЕІМС» ЕК» Еқ2,» RINT» RINC, RN, A, Ils LL, NN, MODE, 
IREND»VALJOsVALJ1»VALJ2»VALNO»VALNI »VALN2»JJ»M» NUMBER» R»I»L»V5K» 
2NsYs YDOT» YL» YLDOT, IIIs COMPARO, COMPAR1,POT,VINT»VINC,MOMENTA, 
3DELTAOs»DELTA1»XSECO,XSECL»TOTAL,RA»ROSsAS»CHANGE ,MANY»S»B» FNORM »BE 
4TTER 


EA EINT 

ОО 199 I = 29311 

ШОШ ЕСІ-1) + EINC 

FK(1) = SQRTF(E(I)/«41499E-24) 


ЕК2(1)= FK(I) * FK(I) 

LL=2*MANY + 1 

πι 

EXPONENTIAL WELL 

FMANY = MANY 

RA(1)=ROS/AS*(1e+F MANY ¥CHANGE) 

70 L=1>LL I 
S=(1.+.-2,.9739%RA(L)+1+,3602%RA(L)%*#2)/(1,-3,7136&%RA(L)+3,1177*RA(L)* 
1*2) 

5555-1. 
В-КА(()%А45/(1.-.904%55%.745%55%%2-.542%55%%3%.174%55%%4) 
POT(L)2S*.751541E-23/8**2 

RN(L)=8/3.5412 


RA(L+1) = RA(L) - CHANGE *ROS/AS 

ШЕ ОС - 1) 190531925197 ; 

БЕГЕТ 196 

FORMAT (13Xs7H ENERGY »s5Xs12H WAVE NUMBER,3X,20H WAVE NUMBER SQUARE 
1D) 


ЕШ 195 + (E(1)> FK(1)sFK2¢1) sI#2sITI) 
FORMAT (3E20.4) 

PRINT 171 

FORMAT(34H EFFECTIVE RANGE/SCATTERING LENGTH) 
PRINT 172»(RA(L) Lz15LL) 

FORMAT (E34. 10) 

PRINT 194. 

FORMAT (13Xs7H RADIUS) 

PRINT 193» ІРМ(1)» іс1911) 

FORMAT (E20.4) 

PRINT 1192 

FORMAT (16H POTENTIAL DEPTH) 

PRINT 2192s (POT(J)sJ=1sJJ) 

FORMAT (F1l602) 

PRINT 191 
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ER 121 
190 


189 
188 
187 


186 
185 
184 


ҒОКМАТ( Зан RADIUS AND ENERGY VALUES COMPUTED) 
DO 187 д ІІ 


ARGB = FK(I) * REND 

CALL BESSEL( АКОВ»2»ЭУАЦМЕ»5» O) 
VALJO( I) = VALUE(Q) 

VALJICI) = VALUE(1) 

VALJ2(1) = VALUE(2) 

CALL BESSEL (ARGBs2sVALUEs 630) 
VALNO(T) = VALUE(O) 

VALN1(1) = VALUE(1) 

VALN2(1) = VALUE(2) 


IF (MODE = 1) 18718735189 

PRINT 18851 | 

FORMAT( 1415) 

CONTINUE 

IF (MODE — 1) 18451865186 

PRINT 185 | 

FORMAT( 26H BESSEL FUNCTIONS COMPUTED) 
RETURN 

END 
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YUKAWA POTENTIAL 


SUBROUTINE READY2 

DIMENSION Е (100) „ЕК( 100) »FK2(100) »RN(100) sVALJO(100) »VALJ1(100) 
1»VALJ2(100)» VALNO(100)» VALN1(100)»VALN2(100)» R(1000)» Y(1000)» 
2YDO0T(1000)» COMPARO(100)» COMPAR1(100)» YL(1000)»YLDOT(1000)»POT(1 
300) »DELTAO(100) »DELTA1(100)»XSECO(100)»XSEC1(100)»TOTAL(100)»RA(10 
40) »FNORM(100)»BETTER(10>100) 

DIMENSION VALUE (40) 

COMMON EsEINTsEINCs FKs FK29 RINT, RINC, RNs, A, Il, LL, NN, MODE, 
IRENDSVALJOsSVALJ1sVALJZ2SVALNOsSVALN1 sVALN2sJJsMs NUMBERS RoslslLsVsK,s 
2м,у, YDOT» YL» УГООТ, 111» СОМРАКО, СОМРАК1,РОТ,УТМТ „МТМС , МОМЕМТА » 
3DELTAO»DELTA1 »XSECO»XSEC1»s TOTAL „КА „КО5 „А5 „ СНАМОЕ , МАМУ , 5 „В, ЕМОКМ „ВЕ 


4TTER 
E(1) = EINT 
DO 199 I - 2,11 
Ge) = Е(І-1) + EINC 
ЕК(І) = SORTF(E(1)/.41499E-24) 


199 FK2(I)s FK(I) * FK(I) 
LL=2*MANY + 1 
= | 
FMANY = MANY 
КА(1)-В05/А5%(1.-ЕМАМҮ%СНАМСЕ) 
YUKAWA WELL 
0 170 ісізіі 
$2(1.—-3.5093*RA(L)-2.0354*RA(L) X*2)/(1647441348*RA(L) 2. 7072 *RA(L ) κα 
l2) 
55-5-1. 
В-КА(1.)%А45/(1.-1.2369%55-4-1.092%55%%2-1.127%55%%2-1.005%55%%4-.,224%5 
15%%5) 
POT(L)2S*.147585E-23*2.1196/8**2 
RN(L)28/2.1196 
170 БА(141) = RA(L) — CHANGE *ROS/AS 
IF (MODE - 1) 19051923197 
197 PRINT 196 
196 FORMAT (13X>7H ENERGY>5X»>12H WAVE NUMBER s3Xs20H WAVE NUMBER SQUARE 
1D) 
PRINT 195 »(Е(І)» ЕК(І)ЭЕК2(1)»91-2,11) 
195 FORMAT (3E20.4) 
PRINT 171 
171 FORMAT( 34H EFFECTIVE RANGE/SCATTERING LENGTH) 
ENT 172,!RA(CL) LS15LL) 
172 FORMAT (E34. 10) 
PRINT 194 
194 FORMAT (13X,7H RADIUS) 
PRINT 193» (RN{tL), L=l;»LL) 
193 FORMAT (E20.4) 
PRINT 1192 
1192 FORMAT (16H POTENTIAL DEPTH) 
BEN NM 21925 (POT(CJ) 923505993) 
2192 FORMAT(F16e2) 
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122 
I9] 
190 


89 
188 
187 


186 
155 
184 


Ро рупа 
FORMAT( 34H RADIUS AND ENERGY VALUES COMPUTED) 
DO 187 1-2,11 


ARGB z FK(I) * REND 

CALL BESSEL( ARGB,25,VALUE»55, O) 
VALJO(I) = VALUE(O) 

VALJ1(1) = VALUE(1) 

VALJ2(I) = VALUE(2) 


CALL BESSEL (ARGB»2sVALUE» 690) 
VALNO(I) = VALUE(O) 

МАЦМ (Г) VALUE(1) 

VALN2(1) VALUE(2) 

IF(MODE - 1) 187,187,189 

PRINT 188; 1 

FORMAT( 1415) 

CONTINUE 

IF (MODE - 1) 18451865186 

PRINT 185 

FORMAT( 26H BESSEL FUNCTIONS COMPUTED) 
RETURN 

END 
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199 


170 


197 
196 





171 
172 
194 
193 

1192 


Е192 
(92 


195 


SQUARE WELL POTENTIAL 


SUBROUTINE READY2 

DIMENSION E(100)sFK( 100) ›»ЕК2 (100) »ь ВМ( 100) » УА ЈО (100) »VALJ1(100) 

WO VALJ2(100)» VALNO(100)» VALN1(100)»VALN2(100)> R(1000)» Y(10007s 
2YDOT(1000)» COMPARO(100)» COMPAR1(100)» YL(1000)»YLDOT(1000)»POT (1 
300) sDELTAO( 100) sDELTA1 (100) »XSECO(100) »XSEC1(100)  TOTAL(C100) 5RA(10 
40) „ЕМОВМ (100) „ВЕТТЕК (10»100) 

DIMENSION VALUE (40) 

COMMON E,EINTsEINCs ЕК» ЕК2, ВІМТ, КІМС» RN, A» II» LL» NN, MODE; 
| ВЕМО, МА: .)О» МАГ.) 1 + МАЦ Ј2 » VALNO» VALN1 »VALN2sJJ9Ms NUMBER» RyI+sl V K, 
2М,Ү» ҮбОТ» ҮШ» ҮШрОТ» ІІІ» СОМРАВО» СОМРАК1»РОТ»ЭУІМТ»УІМС»МОМЕМТА» 
ЗОЕСТАО „ОЕ|ТА1 , ХЗЕСО , ХЗЕС1 „ ТОТА | „БА » ROS »А5 , СНАМОЕ , МАМУ» 5 В, РМОКМ „ВЕ 


4TTER 
END = EINT 
DO 199 I = 2>11 
КОП E(I-1) + EINC 
ЕК(І) - SORTF(E(1)/.41499E-24) 
ЕК2(1) = ЕК(1) * FK(I) 


LL=2*MANY + 1 

JJ=LL 

SQUARE WELL 

ЕМАМУ + МАМУ 

CHANGE SIGN TO - FOR NORMAL OPERATION 

КА( 1) -КОЗ/ИАЗХ (1. +FMANY#CHANGE) 

O ITO L=1>»LL 

5-(1.-1.574%ВҺА(1)-.35З3%КА(1.)%%2)/(1.-2. 384*RA(L)41. ЗОЖКА( |.) #2) 
55-5-1. 

В-КА(() ХА5/(1»- ,5%55%.486%55%%2-.1%16%55%%3%.279%555 Ж ага ME 
End - S*.102276E-23/5**2 


RN(L) =B 

RA(L+1) = RA(L) — CHANGE *ROS/AS 
IF(MODE — 1) 190,192,197 

PRINT 196 


FORMAT (13X>7H ENERGY>»5X>12H WAVE NUMBER»3X»20H WAVE NUMBER SQUARE 
1D) 

PRINT 195 >(E(1)> FK(I)>FK2(1)152,11) 

FORMAT (3E20.4) 

PRINT 171 

ЕОВМАТ ( Зан EFFECTIVE RANGE/SCATTERING LENGTH) 
PRINT 1725 (RA(L) Le1sLL) 

FORMAT ( E34. 10) 

PRINT 194 

FORMAT (13X»7H RADIUS) 

PRINT 193» (RN(L)+s L=1sLL) 

FORMAT (E20.4) 

PRINT 1192 | 

FORMAT (16H POTENTIAL DEPTH) 

PRINT 21929 (POT(J) sJ=l JJ) 

FORMAT (F16.2) 

PRINT 191 я 
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191 FORMAT( 34H RADIUS AND ENERGY VALUES COMPUTED) 
190 00 187 1-2,11! 

ARGB = FK(I) * REND | 

CALL BESSEL( ARGB$25»VALUEs»5, O) 

VALJO(I) = VALUE(O) 

VALIL Cle ae VALUE (ey 


VALJ2(1) = VALUE(2) 
‚ CALL BESSEL (ARGB»2»VALUE» 6,0) 
1 VALNO(I) = VALUE(O) 

VALN1(1) = VALUE(1) 

VALN2(I) = VALUE(2) 


IF (MODE = 1) 187,187,189 
189 PRINT 188,1 
188 FORMAT (1415) 
287 CONTINUE 
IF (MODE = 1) 18451865186 
186 PRINT 185 
185 FORMAT( 26H BESSEL FUNCTIONS COMPUTED) 
184 RETURN 
END 
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GAUSSIAN POTENTIAL 


SUBROUTINE DERIV (AGCsYCsDY) 

BEIMENSION E(100)>FK( 100) »FK2 (I 00995 &N{ 1003, VALJO (780) „VAESITEET OD 
1»VALJ2(1C00)s VALNO(100)» VALN1(100)»VALN2(100)» R(1000)» Y(1000)» 
2YDOT(1000)» COMPARO(100)» COMPAR1(100)» YL(1000)9YLDOT(1000)»POT(1 
300) „ОЕ! ТАО (100) „ОЕГТАТ (100) »ХЖЗЕСО (100) + ХЅЕСІ1 (100) „ТОТАС (100) „КА (10 
40) „ЕМОКМ (100) »ВЕТТЕК (10»100) 

РИЕЕМЕӘТОМ ҮСІЗО)» ОҮ(30) 

COMMON EsEINTsEINCs FKs FK29 RINTs RINCs RNs Aso II, LL, NN, MODE, 
1REND»VALJOsVALJ1»VALJ2sVALNO»VALNIsVALN2,JJ,M, NUMBER» RolIoloVoKs 
2NsYs YDOTs YLs YLDOTs III» COMPARO, COMPAR1,POT,VINT»VINC,MOMENTA» 
3DELTAO »sDELTA1 »XSECOsXSEC1s TOTAL »sRA sROS sAS » CHANGE » МАМҮ› 5, В. ЕМОВМ „BE 
4TTER | 

DY(1) = 

DY(2) = 

RETURN 

END 


Ше(2) 
е ТЕКТЕ + νυ... ο. πο FT πο 


SUBROUTINE DERIVI(XCsYCsDY) 

DIMENSION E(100) »FK( 100) „ЕК (100) „КМ 100) „УА| Ο(100) ΝΑΙ. 1 (100) 

1 »МА. 2100)» МАМО(100)» МА М1(100) „МА М2 100)» К( 1000)» У(1000)» . 
2YDOT(1000)»s СОМРАКО( 100), COMPAR1(100)» YL(1000)»YLDOT(1000)»POT[(1 
300) „ОЕ| ТАО01100) sDELTA1 (100) я ХЅЕСО (100) „Х5ЕСТ (100) » ТОТА! (100) »КА (10 
40) »FNORM( 100) sBETTER(109100) 

DIMENSION YC(300, DY(30) | 

COMMON EsEINTsEINCs ЁК» FK2s RINTs RINCs RNs As» Ils LL» NN, MODE, 
LRENDsVALJO >» VALJ1 » VAL J2 sVALNO s VALN1 sVALN2sJJ9Ms NUMBER» RsI»sL»V>»K; 
2N»Y» YDOT» У|» YLDOTs ІІІ» СОМРАКО» СОМРАВ1»РОТ»УІМТ»УІМС»МОМЕМТА» 
“3DELTAOsDELTA1sXSECOsXSEC1sTOTAL sRA»sROS 5 А5 » СНАМСЕ , МАМУ »5 „В, МОКМ „ВЕ 


4TTER 
ШІСІ) --Үс(2) 
ШЕ! = - (ЕКО (К) + М "““"ЕХРЕ(-(ХС/А)%%2) = 2./ХС%х%2)ХхҮС(1) 
RETURN | Ä : 
END : 
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EXP ON EMA ΟΤΕ... 


SUBROUTINE DERIV ARE Yo DZ 

DIMENSION E(100)»FK(100)»FK2(100)»RN(100) »VALJO(100) ΝΑΙ 11 (100) 
1 «УА1./2(100)» VALNO(100)» VALN1(100) VALN2(100)» R(1000)» Y(1000)» 
2YDOT(1000)» COMPARO(100)» COMPAR1(100)» YL(1000)»YLDOT(1000)»POT(1 
200) „ВЕСТАО( 100) »DELTAL (IE 1 0 DEE ы 
.40) sFNORM( 100) » ВЕТТЕК (10,100) 

DIMENSION ҮС, DY SON 

COMMON EsEINTsEINCs FKs FK25 RINT» RINC», ЕМ. А» ІІ» LL, NN, MODE, 
LRENDs VALJO »sVALJ1 » VALJ2 sVALNO sVALNI1 »VALN2sJJ9Ms МОМВЕК» Б»1>1>У»К» 
2NsYs YDOT» YL» YLDOT» IIIs COMPAROs COMPARL,POT,VINT,VINC,MOMENTA» 
ЗОЕ| ТАО „ОЕ! ТА | „Х5ЕСО , ХЗЕС1 »ТОТА, „КА „КО5, А5 „ СНАМОЕ, МАМУ »5 „В, ЕМОКМ „ВЕ 
4TTER 
DY(1) = 
DY(2) = 
RETURN 
END 


YC (2) 
= (FK2(K) + V * EXPF (- XC/A)) * YC(1) 


SUBROUTINE DERIVI(XC»YC»DY) 
DIMENSION E(100)»FK(100)>»FK2(100)»RN(100)»VALJO (100) »VALJ1(100) 


1»VALJI2(100)» VALNO(100)» VALN1(100) ,VALN2(100)» R(1000) Y(1000)» 


- 2YDOT(1000), COMPARO(100)» COMPAR1(100)» YL(1000)>YLDOT(1000)»POT(1 


300)»DELTAO(100) »sDELTA1 (100) »XSECO( 100) »XSEC1 (100) »TOTAL(100) »sRA(10 
40) „ЕМОКМ (100) „ВЕТТЕВ (10»100) | 
DIMENSION YC(30), DY(30) | 

ΟΟΜΜΟΝ Ε»ΕΙΝΤ»ΕΙ͂ΝΟ» ЕК, Ек?, RINT» RINC, RNs As II» LL, NN, MODE, 
1REND»VALJOsVALJ1>VALJ2»VALNOsVALNI»VALN2sJJsMs NUMBER» R»Is5L5V5Ks» 
2N»Y» YDOT» YL» YLDOT» III» СОМРАКО » СОМРАВ І, РОТ »МТМТ , МУТМС » МОМЕМТА » 
3DELTAOsDELTA1 »XSECOsXSEC1s TOTAL sRA » ROS „А5 , СНАМОЕ , МАМУ »5 „В, ЕМОКМ „ВЕ 


4TTER 
DY(1) = YC(2) 
D2) = = (FK2(K) + V = ЕХРЕСЕ ХС/А) -2. /ХС%%Ұ2)ЖҮС(1) 
RETURN 
END 





YUKAWA POTENTIAL 


SUBROUTINE DERIV (XCsYCsDY) 
BIMENSION E(100)sFK(100) ЕК2 100) RN( 100) »>VALJO(100) sVALJ1 (100) 


1»VALJ2t100)» МАГМО (100) » МАСМ1 (100) „МАСМ2 (100) з R(1000)» Y(1000)» 
2У700ОТ (1000), СОМРАКО (100), СОМРАК (100) » YL(1000) sYLDOT(1000) „РОТ 1 
ESCOSDELTAO(100)»DELTAI(100j),XSECO(100) ,XSECI(C(IDO)S TOTAL CTOO) RACIO 
4&0) „ЕМОКМ (100) „ВЕТТЕК (10»100) 

EDSGNSION YC(30), DY(30) 

COMMON ЕЭЕІКТ»ЕІМС» ЕК» ЕК2, RINT» RINC, RN», A» II» LL», NN, MODE; 
1REND»VALJOSVALJ1sVALJ25»VALNOSVALN15»VALN2» JJ.SM», NUMBER» R»Is5L5,V5,K» 
2Ν9Υ» YDOT» YL» YLDOT, III» COMPARO, COMPAR1,POT,VINT»VINC,MOMENTA; 
3DELTAO,sDELTA1»XSECO,XSEC1»TOTAL»RA»ROS»AS»CHANGE »MANY»S»B» FNORM » BE 


4TTER 5 
Em) = YC(2) қ 
 DY(2) = — (ЕКО(К) + у # (ЕХРЕ(-ХС/А))/ (ХС/А))  УС(1) 
RETURN 
END 


SUBROUTINE DERIV1(XC | YC DY) 

TS STON E(100),FK(100)yFK2(100),RN(100) yVALJ0(—°(100) »sVALJ1(100) 

ο "δι 21100). VALNO(100)» VALN1(100),VALN2(100)» R(1000), Y(1000)» 
2YDOT(1000), COMPARO(100)» COMPAR1(100), YL(1000),YLDOT(1000),POT(1 
LOSE AT e ESO 11007) SECT LADO 1 TOTAL Ы 
40) „ЕХОКМ (100) „ВЕТТЕК (10»100) 

DIMENSION YC(30)» ОУ (30) 

COMMON ЕЭЕІМТЭЕІМКС» ЕК» ЕК2, БІМТ» БІКС» БМ» А» ІІ» LL», NN, MODE, 
1REND»VALJOsVALJ1»VALJ25»5VALNO»VALNI15»VALN2»JJ»M» NUMBER» R»I»5»L»5V5K» 
2NsYs YDOT» YL» YLDOT» III» COMPAROs COMPARI,POT,VINT»VINC,MOMENTA, 
3DELTAO,DELTAI»XSECO»XSEC1»TOTAL»RA»ROS»AS»CHANGE »MANY»S»B» FNORM »BE 


4TTER 

DY(1) = YC(2) 
ШЕШЕ - (FK2(K) + \/ % (ЕХРЕ(-ХС/А)) /(ХС/А) -2./ХС%%2) % YC(1) 
RETURN 
END 


E 





SQUARE WELL POTENTIAL 


SUBROUTINE DERTYV (X€ + rC oD O) 

DIMENSION E(100)>FK1100)»FK21100)»RN(100)»VALJO(100) »VALJ1(100) 
1»VALJ2(100)» VALNO(100)» VALN1(100)»VALN2(100)» R(1000)> Y(1000)» 
2YDOT(1000)» COMPARO(100)» COMPAR1(100)» YL(1000)»YLDOT (1000) »POT (1 
300) „ОЕГТАО (100) „ОЕ| ТАТ (100) »ХУЕСО (100) » ХЅЕС1 (100) »ТОТАГ (100) »КА (10 


40) „ЕМОКМ (100) „ВЕТТЕК (10»100) 

DIMENSION YC(30), DY(30) 

COMMON E»*EINTsSEINC$, FK, FK25 RINT» RINC, RN» A» II» LL» NN», MODE, 
IREND»VALJO»VALJ1»VALJ25»VALNOsSVALNI15»VALN25»JJ»M, NUMBER» R»I»L»V5»K, 
2NsYs YDOTs YLs YLDOTs IIIs COMPAROs COMPAR1 »POT  sVINT  sVINC sMOMENTA 9 
3DELTAO »sDELTA1 »XSECOsXSECL1lsTOTAL»sRA »ROS 9 AS 5 CHANGE » MANY S989 F NORM „ВЕ 
4TTER 


Dri l) =YC(2) 
IF(XC -Α) 1»1»2 


ES(2) - -(FR2K) t V) * YC (1) 
"RETURN 

2) ==(FK2(K))%* YC(1) 
RETURN 

END 


SUBROUTINE DERIV1(XC YC DY) 

DIMENSION Е(100) »ЕКС 100) »9FK2(100) sRN( 100) »VALJO(100) »VALJ1(100) 
1»VALJ2(100)5» VALNO(100)» VALN1(100)5,»VALN2(100)» R(1000), Y(1000)» 
2YDOT(1000)» COMPARO(100)» COMPARI1(100), YL(1000)sYLDOT(1000) sPOT(1 
300) »ОЕ: ТАО (100) „ОЕ: ТА (100) „ХЗЕСО (100) »Х5ЕСТ (100) „ТОТАГ (100) »КА (10 
40) „ЕМОКМ (100) ЭВЕТТЕК(10>100) | 

EIMENSION YC(30)5 DY(30) 

COMMON ES,EINTsSEINC, FK,» FK2» RINT, RINC, RN, As» II» LL» NN», MODE, 
1REND»VALJO»VALJ1»VALJ2»VALNO»VALNI»VALN2»JJ»M» NUMBER» RsI»L»V5K» 
2N»Ys YDOT, YL» YLDOT, III» COMPARO, COMPARI,POT,VINT»VINC,MOMENTA;> 
3DELTAOsDELTAI»XSECOs,XSECI»TOTALsRAsROSsSASsSCHANGES MANY 3 5, B» FNORM BE 


4TTER 
DY(1) =YC(2) 
MAC =A) 1,1»2 


1 0М(2) + -(ЕК? (К) + V - 2./ΧςΧΧ2)ΧΥς(1) - 
RETURN | 

2 DY(2) s -(FK2(K) - 2./ XC**2) *YC(1) 
RETURN 
END 
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mee FOLLOWING DEFINITIONS FOR B MUST BE SUBSTITUTED IN READY2 FOR УАКҮТИС 
Meee RING LENGTHe THE DEFINITION PROVIDED IN READY2 IN THE LISTING IS FOR 


ARIATIONS IN EFFECTIVE RANGE. 


YUKAWA WELL 
ο ο ο ООШ s 
SQUARE WELL 
В-Қ05/(1.-.5%(5-1.)4%.486%(5-1.)%%2-.476%(5-1.)%%3-.378%(5-1.)%%д 
ЕЕ. 1208%(5-1.)%%5) 
EXPONENTIAL WELL 
„ВеНО5/(1.-.904#55+.745#66##2-„543#5би#3+.174#59##4) 


GAUSSIAN WELL 


В-К05/(1.-.641%55%.568%55%%2-.721%55%%3%1.172%55%%4-1.015%55%%5) 
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These results are for a yukawa potential with Г аа figed at 
2 


2.46 Fm and a. varying in | 16% deviation from a¿* = 23.677 FM. 
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